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Abstract. We propose an approach to the aggregation of risks which is based on estima- 
tion of simple quantities (such as covariances) associated to a vector of dependent random 
variables, and which avoids the use of parametric families of copulae. Our main result 
demonstrates that the method leads to bounds on the worst case Value at Risk for a sum of 
dependent random variables. Its proof applies duality theory for infinite dimensional linear 
programs. 



1. Introduction 

Aggregation of risks is a key issue for risk management and regulation in the financial 
sector. Consider an institution which may incur losses Xi, . . . ,X„ in each of its n differ- 
ent business divisions. To assess its overall exposure to risk, knowledge of the distribution 
of the sum S = Xi + . . . + Xn is vital. Suppose that the distributions of the individual 
Xi are reasonably well-known. This does not imply, however, that there are enough data 
available to estimate the distribution of S, let alone the joint distribution of the vector 
X = {Xi, . . . , Xn). As the process of implementing the Basel II and Solvency II accords has 
amply demonstrated, the data on which the knowledge of the marginal distributions is based 
are often not suitable for reconstructing the joint distribution: They may have been collected 
at different times, asynchronously, or subject to wildly differing standards of exactness or 
documentation (see, e.g., [IS])- Nonetheless, it is often the case that at least some data from 
the joint distribution are available — perhaps not enough to obtain a reasonably accurate 
picture of the tail behavior of S, but sufficient to estimate quantities such as the covariances 
of the Xi. 

Suppose now, for simplicity, that one is ultimately interested in the Value at Risk of S 
(at some level a). If there were no information whatsoever available about the dependence 
structure of X, one would have to resort to a worst case analysis in the sense that one would 
let the joint distribution of X run through all n-variate distributions with the appropriate 
marginals. In the situation at hand, with at least some aspects of the dependence structure 
known, it is natural to take only those joint distributions into account that are compatible 
with the additional information. 

Arguably, this procedure makes sense even in the case when one has enough data from 
the joint distribution to feel comfortable with estimating the distributions of S and even 
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of X. In this situation, the worst case analysis subject to a constraint might be viewed as 
a robust version of VaR estimation (robust in the sense of taking model uncertainty into 
account, see [5J) — even more so as the constraints we are concerned with are stated us- 
ing quantities that are amenable to robust estimation (in the statistician's sense of the word). 

The present paper treats the "worst case value at risk" of a sum of two random variables 
whose covariance is assumed to be known. We rephrase the worst case analysis as an infinite 
dimensional linear program and establish a handy dual version, from which it is easy to 
extract bounds for the worst case value at risk. Making the tools of duality theory bear 
upon the present setup constitutes the technical core of what follows. 

Before providing, in Section [2], an exact statement of the main result, let us try to place 
our approach within the context of risk management literature and practice. The last decade 
has seen an upsurge of interest in the concept of copulae. Parametric families of copulae 
have been widely used as a basis for estimation of various quantities related to aggregate risk 
(see, e.g., pj and the references therein), even as these procedures have also met fundamental 
criticism (see [IDj)- By contrast, the methods which are advertised in the present paper 
avoid any use of parametric families - as well as of inversion of marginals - and are based on 
estimation of scalar or finite-dimensional quantities, not of real valued functions on a unit 
cube. 

Among the literature that is based on the copula concept, it is a paper of Embrechts, 
Hoing, and Juri [1] that is closest in spirit to the present approach. The authors use order 
relations for copulae to exploit partial information on the dependence structure of X for 
worst case bounds on VaRa(5'). Specifically, if C denotes the copula of X and if there are 
copulae Co, Ci such that Cq < C and < Cf (C^ denoting the dual copula), then the 
authors find upper and lower bounds for VaRa(S'). The copulae Cq, Ci may be chosen to 
encode properties of X such as positive orthant dependence. In contrast to that, the present 
approach is more empirical in spirit, in that the partial information we use can be easily 
estimated from the data available. 



Suppose that X and Y represent losses associated with different positions, and that a risk 
manager wishes to calculate VaRa(5'), the Value at Risk at the level a, oi S = X + Y , that 
is, the infimum over all z G M such that 



where /i denotes the joint distribution of the vector (X, Y). Suppose now that nothing about 
yU is known, except for that its marginals are P, the distribution of X, and Q, the distribution 
of F, respectively. Then it seems natural to replace VaRa(5') with WVaRa(P, Q), which we 
define as 



2. Worst case analysis based on marginals and covariance 





inf 

fieM{P,Q) 




with 



M{P, Q) := {fi E M\R'') : fi o pi-' = P, /i o pr2 ^ = Q} . 
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Here pr^^ and pig are the projections on the components of M x M, and A^^(R^) denotes the 
set of probabihty measures on M^. Note that WVaRa(P, Q) > VaRa(S'). It is thus natural to 
regard WVaRa{P,Q) as a "worst case Value at Risk", taking care of the model uncertainty 
inherent in estimating quantiles of the distribution of a sum of dependent random variables 
whose joint distribution is not (completely) known. 

Now suppose that some (though not necessarily exceedingly copious) data from the joint 
distribution of {X,Y) are available. These can be used to estimate, e.g., the covariance 

xy fi{dx, dy) - J xP{dx) j yQ{dy), 

which typically can be estimated in a robust way on the basis of rather few data - in stark 
contrast to the wealth of data that would be necessary for a reliable estimate of the copula 
of {X,Y). Incorporating this extra information into the worst case analysis will lead to a 
restricted worst case Value at Risk that will in general be smaller than WVaRa{P,Q). To 
be specific, set 

WVaRa(P, Q, k) := int I z e R : inf / / lf„oo z)ix + y) u(dx, dy) > a> , 

{ ,ieM{P,Q,k) J J ' J 

where 

M{P,Q,k):=\^ieM{P,Q): [[ xy fi{dx,dy) = k\. 



Note that the obvious choice for A; is A; = px,Y + epCg, where px,y is an estimate of the 
covariance of X and Y, and 



ep := J xP{dx) and eq := j yQ{dy) 



are assumed to be known. 



The crucial problem, which will be treated in what follows, is how to evaluate an infimum 

^l&M{P,Q,k) J J 

where i{x,y) is a lower semicontinuous bounded function such as l(_oo,2)(a; + y). (This level 
of generality should provide some leeway for discussions of more general aggregate financial 
positions (see [HJ p. 248]) and risk measures.) In the absence of extra information about the 
joint distribution, i.e., if the infimum in ([1]) is taken over Ai{P, Q) rather than Ai{P, Q, k), 
the problem is but an instance of the classical (Monge-Kantorovich) mass transport prob- 
lem: A substance whose initial spacial distribution is given by P, is to be shipped to a final 
spacial distribution given by Q, according to a transportation plan which can be encoded 
by an element of Ai{P,Q) and which is supposed to minimize the overall cost if transport 
from X toy comes at the price i{x, y). Monographic treatments of mass transport are due to 
Rachev and Riischendorf [Ill|12] and Villani [15]. The special case of i{x, y) = l(_oo,z)(a; + y) 
has been studied in Makarov [8], Riischendorf [H], and Frank, Nelsen, and Schweizer |6]. 
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The modern functional-analytic approach to mass transport, which is due to Kantorovich, 
is based on the following duality: 



// ^{x,y)l^{dx,dy) = sup / f{x)P{dx)+ / g{y)Q{dy) ] , (2) 

where 

Af{^) := {(/, g) G a(M) X a(M) : f{x) + g{y) < y) Vx, y G M} , 

C6(]R) denoting the set of bounded continuous functions on M. Note that by plugging suitable 
test functions into the right-hand side, one obtains lower bounds on the left-hand side. Re- 
turning to WVaRa(P, Q), one observes that replacing for each z and i{x, y) = l(_oo,z)(x + y) 
the left-hand side of ([2]) by a lower bound will lead to an upper bound on WVaRa{P,Q), 
and this is what matters in risk management. (Of course, a lower bound on WVaRa(P, Q) 
can be obtained by plugging, for each z, test measures into the left-hand side.) 

Given the usefulness of ([2]), one would like to have an analogous statement for the case 
that the infimum is taken over Ai{P, Q, k) rather than A4{P, Q). This is the content of the 
main result of this paper, which will be proven in Section [S] 

Theorem 2.1. Assume that i is bounded and lower semicontinuous and that P,Q E A4^(R) 
satisfy 

J x^ P{dx) < oo, j y"^ Qidy) < oo. 

Then for all k eW with 



inf xy fi{dx,dy) < k < sup xy ^{dx,dy), (3) 

one has that 

inf / / i{x,y) fi{dx,dy)= sup I / f{x)P{dx) + / g{y)Q{dy) + ak 1, (4) 

t^&M[P,Q,k) J J {f,g,a)eU{e) 

where 

Uie) := {(/, g, a) G S(R) x ^(M) x M : f{x) + giy) + axy < i{x, y) Wx, y eR} 

and 

S{R) := {/ G C(M) : f{x)/{x^ V 1) is bounded} . 
Moreover, a minimizing measure E A4{P, Q, k) for the left-hand side of (jl]) exists. 



Remark 2.2. 

(1) We note that if the supports of the measures P and Q are bounded, then in the 
definition of U{i) one can replace S{R) x S{R) by Cb(supp(P)) x C;,(supp(Q)). 

(2) It is not clear whether there exist maximizing vectors (/, g, a) for the maximization 
problem on the right-hand side of (jl]). In order to obtain such kind of result one 
usually needs apriori bounds on /, g and a. If a is bounded, apriori bounds can be 
obtained by replacing / and g by suitable new functions /o and go (see, e.g., 
However, in our situation no apriori bound for a is available. 
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(3) To obtain explicit values for the bounds on in ([3]), note that they are solutions to 
classical transportation problems with cost function c : — M; (x, y) ^ xy. The 
cost function c is twice continuously different iable and the mixed partial derivatives 
Cxy are positive. Therefore there exists an optimal transportation plan for those opti- 
mization problems; i.e.: there exist kernels T~^{x, dy) and T~{x, dy) with the property 
that the measures fi~^{dx,dy) := T'^{x,dy)P{dx) and ^~{dx,dy) := T~{x,dy)P{dx) 
satisfy 

inf xy ^{dx,dy) = xy fi'{dx,dy) 

ti&M{P,Q)JJ J 

and 

sup xy fi{dx,dy) = xy fi^{dx,dy). 

li&M{P,Q)JJ J 
In order to provide a more explicit representation for those kernels in terms of the dis- 
tribution functions F and G of the measures P resp. Q, we introduce the generalized 
inverse of a non- decreasing right-continuous function if : M ^ M as 

H-{i) ■= mi{x e M+ : F{x) > 

the left-limit of H at the point x G M as H{x—) := lim^j^o H{x — e) and the jump 
of at a; G M as AH{x) := H{x) — H{x—). Further, we use the notation 6u{dv) 
to denote the Dirac measure on M, which assigns a unit point mass to a single point 
M G M. Moreover, for a Borel-measurable set A C M with positive Lebesgue measure 
we denote by Uji{dy) the uniform distribution on A. We then have 

r MA ._ / SG-{F{x)){dy) if AF{x) = 

L+[x,ay) . I f/[G-(^(._)),G-(F(.))](cil/) if AF{x) + 

and 

r d„\ ■= [ ^G-{i-F{.)){dy) if AF{x) = 

i.[x,ay). \ u^^-^^^F(.~)),G~il-Fi.))]{dy) if AF(x) ^ " 

This explicit representation follows from the fact that the supports of the optimal 
measures satisfy strong constraints (see [I], p.96]), which restricts the set of possible 
transportation plans to marginal transformations, which transform the measure P to 
the measure Q. 

In the case that X or y has a heavy-tailed distribution. Theorem 12.11 is not applicable as 
stated. In this situation, one may replace xy by {xy V —R) A R for a suitable R, a natural 
choice from the point of view of robust statistics. Then the following variant of Theorem 12. II 
applies, whose proof is actually a simplified version of the one given in Section [3] for Theorem 
12. II and will therefore be omitted. For n G C;,(]R^) write 

M^{P, Q, k) := |/i G M{P, Q) : jj k{x, y) fi{dx, dy) = k 

Theorem 2.3. Assume that i is bounded and lower semicontinuous and n G C;,(M^). Then 
for all k with 

inf / / K{x,y) fi{dx,dy) <k< sup / / K,{x,y) fi{dx,dy)^ 

ti&M{P,Q)JJ i,eM{P,Q)JJ 
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one has that 

inf / / £{x,y) fi{dx,dy)= sup I / f{x)P{dx)+ I g{y)Q{dy) + ak \ , (5) 
ti€M4P,Q,k) J J {f,g,a)eu{e) \J J J 

where 

W := {(/, g, a) E CbiM?) x M : f{x) + g{y) + y) < i{x, y) Vx, y G M} . 

Moreover, a minimizing measure /i G ^A^^{P, Q, k) for the left-hand side of ([5]) exists. 

3. Proof of Theorem 12.11 

In this section we will apply the theory of linear programming in infinite dimensional 
spaces (see, e.g., Anderson and Nash [Ij) and will thus have to rephrase our problem in 
the framework of topological vector spaces. Denote by 7^b(M^) the space of signed Radon 
measures with finite total variation norm on M^. TZh{M.'^) contains the cone 7^+(M^) of finite 
nonnegative Radon measures. By the Hahn decomposition theorem, every signed measure 
/i G 7^fe(]R^) can be decomposed into /x = /i+ — where //+ and /i~ are elements of 7l+(M?), 
and one can define the total variation measure := + /i~ (see, e.g., [2]). Consider the 
vector spaces 



X := <j /i G 7^^,(M^) : / / + r d\l^\ < oo 
and 

Y:={hE C(R2) : h{x,y)/iix' + y') V 1) G ^(M^)}. 
Note that Y contains the function (x, y) ^— xy. The pairing 

(., .) : X X Y ^ M; (yU, /i) ^^ (fi, h) := // h{x, y)fi{dx, dy) 



puts X and Y into duality. In X one has the cone /C := X fl 7^+(M^). We endow X with 
the cr(X, Y) -topology, which is the coarsest topology such that for all /i G Y the functionals 
fl ^ {fx,h) are continuous. Observe that Cfe(M^) C Y, so cr(X, Y)- convergence implies weak 
convergence in the usual measure-theoretic sense. (In what follows, "weak concergence" 
(written as — — ->■) will always be understood in this sense. Otherwise, the relevant a-topology 
will be made explicit.) We collect this and related useful properties in the following 

Lemma 3.1. We write Vk ■= {{x'^ + y'^) V l)/i„, u := ((x^ + y'^) V l)/i. 
(i) 

tr(X,Y) w 



m 



(iv) 



/in > /i =^ [l^n) IS tight. 



{vn) is tight =^ (/in) is tight. 
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We define the vector spaces 



and 



j d\u\ < oo| X |z/ G TZbiR) : j d\v\ < oo| x 



/ G cm : f{x)/{x^ V 1) G x j^? G C(M) : g{y) / {y^ V 1) G CWj x 

The vector spaces V and W form a dual pair with pairing 

(.,.): VxW^M; ((/ii, /is, a), (/, ^7, J f{x)fii{dx) + j g{y)i,2{dy) + ap. 

Define the constraint operator 

v4:X->V; fi ^-^ /io pr^^^ o pr2 \ xy fi{dx, dy) 



The problem to evaluate the left-hand side of (jl]) now can be restated in linear program- 
ming language as follows: 



minimize: JJ £{x,y)fi{dx, dy) (6) 

subject to : Afi = (P, Q, k) 
fi e IC. 

We will refer to ([6]) as the primal problem. To see its equivalence to the original problem, 
note that square integrability of the marginals implies G /C, and that finite Borel measures 
on are Radon measures (see, e.g., |2]). Measures G /C which satisfy the constraints in 
the primal problem are called feasible. Note that such /i must be probability measures. If 
feasible solutions to the primal problem exist, it is called consistent. The minimal value of 
the integral in is called the value of the primal problem in optimization theory. In what 
follows, we will use the more suggestive terminology of "optimal value" . 

In order to state the corresponding dual problem, we have to compute the adjoint operator 
A* -.W^Y oiA. We have 



{Afi,{f,g,a)) = { ifiopT^^,fiopT^\ xy fi{dx,dy)] ,{f,g,a) 

J{x) + g{y) + axy^fi{dx, dy). 

Thus we see that 

A* : m 



Y; (/, g, a) ^ (^(x, y) f{x) + g{y) + axy^ . 
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The dual problem is now given by 

maximize : J f{x)P{dx) + j g{y)Q{dy) + ak (7) 

subject to : A*{f,g,a.) < £ 
{f,g,a)eW. 

Triplets {f,g,a) which satisfy the constraints of the dual problem ([71) are called feasible. 
The dual problem is called consistent, if feasible triplets exist. 



The central concern of linear programming in infinite dimensional spaces is to give sufficient 
conditions for the optimal values of the primal problem ([6]) and the dual problem ([7]) to 
coincide. The weak duality theorem (see [U Thm.3.1]) states: 

Proposition 3.2. // the primal problem ^ and the dual problem ^ are consistent, then 
the optimal value of the primal problem is larger than the optimal value of the dual problem. 

The difference between both values is called the duality gap. So we have to prove that 
there is no duality gap in the present situation. The following proposition collects some 
essential facts from infinite dimensional linear programming. 

Proposition 3.3. Assume that the set 

H := +r) : G /C, r > o| 

is closed m V x M. Then one has: 

(a) // the primal problem ^ is consistent with finite optimal value, then there is no 
duality gap. 

(b) If the primal problem ^ has finite optimal value, then it is solvable, i.e., a minimizing 
/i G /C exists. 

Proof. These are Theorems 3.9 and 3.22 in Anderson and Nash [1], plus the observation that 
the proofs of these theorems as well as the proof of [H Thm.3.3] do not use the continuity 
of the functional ^— (/x, £). □ 

Lemmata 13.41 and 13.51 will help us prove that in the present set-up the set H is closed. 

Lemma 3.4. For a bounded lower semicontinuous function : — > M, weak convergence 
in 7^+(M^) of fin towards /i* implies 

liminf // h{x,y)fin{dx, dy) > / / h{x,y)fi^:{dx, dy). 



n— >oo 



Proof. We may assume without loss of generality that h is non-negative. We define the 
functions 



m 

k=l 



Note that this is actually a finite sum of indicators of open sets. It then follows that 

h--<g*rr.<h. 

m 
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Using the Portmanteau theorem (see, e.g., [3| p. 24]) one can deduce from this that 

hminf / / h{x,y)finidx, dy) > hminf / / g^{x , y) ^n{dx , dy) 
J J J J 

> J J gm{^,y)f^*{dx,dy) 

> // h{x,y)^^,{dx,dy) — 1/m 



Since m G N is arbitrary, the resuh follows. □ 

Lemma 3.5. The sequence of weighted measures ((x^ +1/^) V is tight if, and only if, the 
two sequences (x^ V l)(yUri ° P^i^) (^i^d (y^ V o W2^) ^^^^ ^^^^ tight. 

Proof. It is obvious that tightness of (x^ + y"^) V 1 /i„ implies tightness of (x^ V o pr|f ^) 
and V o pr^^). 

In order to prove the converse we note that for a prescribed e > there exists > 1 with 
the properties 

/x^ V 1 /i„ o Y>T:i^{dx) < t and / V 1 /i„ o prg ^(rfy) < e. 

It follows from this that 

zl /i„ o pr^^([-2;o, ZqY) < e and zl /i„ o prg ^([-^o, ^o]'') < e- 
This implies that 

(x^ + ?/^) V 1 fin{dx,dy) 

{[-zo,zq]x[-zo,zq])'^ 

x^ V 1 /in(cix, dy)+ 2/^ V 1 /i„(cix, dy) 

-20,20]X [-20,^0]'' >^ ^ [-20,20] X [-20,20]" 

+ / / x^ V 1 yu„((ix, dy)+ Z/^ V 1 /i„((ix, dy) 

^ J [-20,2o]'=X[-20,2o] J ^[-20,2o]'=X[-20,2o] 



+ / / X y 1 finidx.dy) + // y V 1 /i„((ix, dy) 

'[-20,2o]'=X [-20,20]= ^ J[-20,2o]'=X[-20,2o]'= 



< z^finiM X [-^0, zof) + I V 1 /x„ o pr2 ^(d?/) 

"20, 20]= 



+ / x^ V 1 /i„ o pr^ ^(rfx) + 2o/i„([-2;o, ZqY x R) 

[-20,2o]<= 

+ / x^ V 1 /i„ opr^^(rfx) + / V 1 /i„opr2'^((i?/) 

^[-20, 20]^= J[-2o,2o]'= 

< 6e. 

The tightness of the weighted sequence ((x^+?/^)Vl)/i„, follows from these considerations. □ 
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The next lemma will be useful to prove the consistency of the primal problem (Q. It 
concerns the following transportation problems: 



minimize : J xy ^{dx, dy) (8) 

subject to : fio prj"^ = P, /i o pr2^^ = Q 

and 

maximize : J xy jj,{dx, dy) (9) 

subject to : fio pr^^ = pi^^ = Q 

fi E K,. 

Lemma 3.6. There exists an optimal solution fi~ for the optimization problem ^) and an 
optimal solution /i"*" for the optimization problem 

Proof. It suffices to indicate the proof for ([8]). We first note that the product measure 
Hq := P Q satisfies /iq o pr]"^ = P and /xq ° P^2^ = Q- Thus the transportation problem ([8]) 
is consistent. Moreover, we have that 

xy dy) 



< jj{x' + y')yi fi{dx, dy) 

< y a;^ V 1 yU o pT^^(dx) + J V 1 yU o pi^^^dy) 

< J x^yi p{dx) + j y^'^'^ Q{dy) < ^■ 

Thus the transportation problem ([8]) has finite optimal value. To prove the existence of an 
optimal solution, one apphes an analogue to Proposition 13.31 (b) (see [H Thm. 3.22]), which 
states that it is sufficient to prove that the set 

I,:^{(.op.,..,op.,-//.„„,,.)):.e/c} 

is closed in (t(V, W). Assume that there exists a sequence /i„ in /C with the property that 



converges with respect to cr(V, W) towards (/ii, fi2,c). The convergence of these sequences in 
the (t(V, W)-topology implies that the sequences (x^Vl)(/i„opr^^) and (?/^Vl)(/i„opr2 ^) are 
tight. By Lemmata l3.5l and l3.ll it follows that the sequences z/^ := ((x^+?/^)Vl)yU„ and fin are 
tight. Hence there exist subsequences and z/^^. converging to fi^:, z/* = ((x^ + y^) V 
respectively. This means that /z^^, cr(X, Y)-converges to fi^. It follows from this that the 
three sequences 

Hn o pr]"\ Hn o pr2"^ and xy Unidx, dy) 
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converge respectively towards 

yU*opr]^\ fx^opr'^ and JJ xy jj,^{dx,dy), 

hence 

/ii = yU^, o prj"\ fl2 = fJ'*°W2^ c = JJ xy fi^{dx,dy). 

Thus the set D is closed and the proof of ([H]) is complete. □ 

Proof of Theorem \2.1[ We first prove that the primal problem ([6]) is consistent if A; G M 
satisfies 



inf xy jji{dx^dy) <k< sup xy jj,{dx,dy). 

tMeMiP,Q)JJ tMeM{P,Q)JJ 

In Lemma [3l6l we saw that the supremum and infimum are attained in /i"*" resp. fi~ . So there 
exists A G [0, 1] with the property that fix := A/i"*" + (1 — A)/i~ G satisfies 

fixopr^^ = P, fix°PT^2^ = Q and xy fix{dx,dy) = k. 



This proves that the set of feasible solutions for the primal problem ([6]) is nonempty, hence 
is consistent. Moreover, since £ is bounded, the objective functional 



is bounded on the set of feasible solutions. This implies that the optimal value must be finite. 

From Proposition 13.31 it follows that in order to prove Theorem 12. 11 it is sufficient to show 
that the set 

:= +r) : /X G /C, r > o| (10) 

is closed in V x R. 

Consider a sequence of measures /i„ G /C and a sequence of real numbers r„ > such that 
the triplets 



/i„opr/, /i„opr2\ J J xy fin{dx,dy) 
converge in cr(V, W) towards a triplet (/ii, /i2, c) G V, and that 

i{x,y)fin{dx,dy) + r„ 



converges to a real number b. We have to find a measure /x* G /C and a real number > 
with the properties 



/i*opr^ = /xi, /i*opr2 = /i2, J J xyfi^{dx,dy) = c and J J £{x,y)fi^,{dx, dy) + = b. 

The convergence of these sequences in the o"(V, W)-topology implies that the sequences 
(x^ V l)(/in o P^r^) and (y^ V l)(/i„ o pr^^) are tight. By Lemmata 13.51 and 13.11 it follows 
that the sequences z/„ := ((x^ + y"^) V l)/i„ and /i„ are tight. Hence there exist subsequences 
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firik and converging to /i*, z/^ = ((x^ + y"^) V respectively. This means that fin^ 

cr(X, Y)-converges to fi^,. It follows from this that the three sequences 



/i^opr/, Unopi^^ and J J xy Hn{dx,dy) 
converge respectively towards 

/i*opr^\ /x^,opr2^ and JJ xy n^{dx,dy), 

hence 

/2i = /i*opr^\ /i2 = /i* o prj"^ and c = JJ xy n^:{dx,dy). 
Moreover, since i is bounded and lower semicontinuous, by Lemma 13.41 one has 
liininf jj i{x,y)iJ,n^{dx, dy) > jj £{x, y) n^{dx, dy) > -oo. 
By further thinning out the subsequence , we may assume without loss of generality that 

hrn jj i{x,y)fin^,{dx,dy) > jj i{x,y)fj.^{dx,dy) > -oo. 
Together with the convergence of the sequence 

i{x,y)nn{dx,dy) + r„, 

this implies that the sequence r^^ converges to a non- negative real number Too- If we define 

:= Too + hm jj i{x,y)nnk{dx, dy) - jj i{x , y) n^{dx , dy) , 

it follows that 

lim ( 1 1 i{x,y)fin{dx,dy) + r„ ) = ^lirn ( / / i{x,y)finSdx,dy) + 

£{x,y)fi:,{dx,dy) + r*. 

Since by definition r* is non-negative, it follows that the set H is closed. This completes the 
proof of Theorem 12. 1[ □ 
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